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Abstract 

A relativistic non-local model for the inclusive photoproduction of 
^ \ T) mesons from complex nuclei is introduced. The model is based on 

the dominance of the Sn(1535) resonance. We compare the results of 
our calculations with the available data on inclusive cross sections for 
k> , the nuclei C, Al and Cu. Assuming the resonance propagates freely in 

!-h ' the nuclear medium, we find that the calculated angular distribution 

and energy dependence of the cross sections reproduce the data in a 
reasonable fashion. The present non-local model allows the inclusion 
of density dependent mass and width in the calculations. Including 
these in the calculations reveals that the presently available data do 
not show clear preference for the inclusion of such modifications of the 
properties of the Sn(1535) in the nuclear medium. 
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1 Introduction 

The question of possible modifications of hadron properties in the nuclear 
medium has attracted considerable interest in recent years. This interest 
stems from the hope that studies of these effects can clarify the interplay 
of the competing degrees of freedom in the nucleus. In particular it may 
be possible to learn about explicit quark degrees of freedom in the nucleus 
through such studies. 

There have been several theoretical as well as experimental investigations 
of these effects, and evidence is gathering that indeed such modifications 
do take place. Several studies have been carried out for the modifications of 
meson properties in heavy ion and other collisions [1, 2, 3, 4, 5]. Photonuclear 
reactions also lend themselves to investigations of this type [6, 3, 7]. In 
particular, a promising area of investigation is the photoproduction of rj 
mesons on nuclei in the energy range corresponding to the second resonance 
region of the nucleon. The fact that one resonance, the Sn(1535), dominates 
the production process at these energies, makes it attractive to look at these 
reactions for clues on changes in its properties in the nuclear medium. 

Three pioneering and complementary measurements of the inclusive eta 
photoproduction on nuclei have been reported in recent years. The earlier 
measurements of Robig-Landau et al [8] at MAMI included cross sections 
on a number of nuclei from threshold up to a photon energy of 800 MeV. 
These cross sections displayed the expected rise on the low energy side of 
the resonance. Few years later measurements of the same reaction were 
performed at KEK [9], extending the energy range to lGeV. This group 
reported that a broad resonance due to the excitation and decay of the S n 
resonance in the nucleus has been observed for the first time. In a recent 
investigation, Kinoshita et al [10] reported on measurements of this reaction 
on C and Cu targets at LNS, for photon energies up to 1.1 GeV. These 
authors indicated the importance of the contributions from other nucleon 
resonances in the second resonance region in addition to the dominant Sn, 
at energies higher than 900 MeV. The authors of references [9, 10] performed 
calculations for their data using an adaptation of the Quantum Molecular 
Dynamics model. Lehr et al [11] performed calculations based on the semi- 
classical BUU transport model. Another approach was that of Maruyama 
and Chiba [12] who carried out relativistic calculations in infinite nuclear 
matter based on an effective Lagrangian and a scaling factor to account for 



the final state interactions of the 77 mesons. 

In the present paper we present an analysis of the data based on a rela- 
tivistic model for the inclusive photoproduction on finite nuclei. The main 
ingredients of model are the use of an effective Lagrangian approach (Ben- 
merrouche et al[13\) and the relativistic mean field approach to nuclear dy- 
namics [14]. The final state interactions of the outgoing 77 meson are taken 
into account. This fully non-local model has been discussed recently for the 
cases of exclusive and incoherent photoproduction [15, 16]. These non-local 
calculations, in contrast to earlier local ones [17], allow us to include possible 
density dependent changes in the mass and width of the resonances in the 
nuclear medium. 

The justification for the claim that the inclusive data from KEK, Mainz 
and LNS may reveal changes in the properties of the Sn resonance requires 
some discussion. It is known that the Sn is the main contributor to the 77 
photoproduction cross section on free nucleons in the 0.6 - lGeV range. There 
are contributions from neighboring resonances and from nucleon Born terms 
as well as meson exchange diagrams; but these are small by comparison. We 
will show in the course of our discussion that this situation carries over to the 
inclusive production on nuclei. The contributions are only at the few percent 
level, a situation which allows us to probe whether the existing data does 
reveal any information regarding the properties of the Sn in the medium. 
As we shall show there are remaining uncertainties regarding other aspects 
of the analysis which are more important compared to the small contributions 
of other diagrams. 

The plan of the paper follows a recent discussion of the non-local rela- 
tivistic model for the exclusive quasifree photoproduction of 77 mesons. The 
reaction amplitude is cast in a general form in section 2 where the method 
of evaluation of the amplitude in its non-local form is outlined. The results 
of the model and the discussion are given in section 3. Conclusions are given 
in section 4. 



2 Formalism 

In the inclusive reaction, a photon is absorbed by the nucleus and as a result 
an 77 meson is produced. This meson is the only detected particle and there 
is no information about the final state of the nuclear system. In the present 



work it is assumed that the inclusive reaction is based on the exclusive re- 
action in which a single nucleon is ejected in the continum in the course of 
the production process [18, 19]. The transition amplitude for the exclusive 
reaction through the Sn resonance, in its full non-local form, is given by [15], 
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where En, uj v and u; 7 are the energies of the final nucleon, 77 meson and 
incident photon, respectively. M is the nucleon mass, Mr and T represent 
the mass and width of resonance. S^jAJb) and (Jj, J B ; Mf, M B \J t , Mj) are 
spectroscopic and Clebsh-Gordon coefficients. Kr and g V NR are the anoma- 
lous magnetic moment of the resonance and its coupling constant with r\ and 
nucleon. i/) s f(y), ^B^y) and (f> v (y) are wave functions describing the final 
nucleon, initial nucleon and rj meson, respectively. ip B {y) is a solution of the 
Dirac equation with appropriate mean field potentials and <p v is a solution 
of the Klein - Gordon equation. As is explained in [19], the wave function of 
the final nucleon (ip s f(y)) is taken to be a plane wave. Calculations without 
inclusion of the final state interaction of r\ meson revealed that the u-channel 
contribution (second term in (1)) is less by two orders of magnitude than 
that of s-channel (first term in (1)). As a result the non-locality and possi- 
ble medium modifications of the former will have a negligible effect on the 
calculated cross section. We thus ignore this very small change and treat the 
u-channel locally as in [19]. The dominant s-channel diagram is treated in 
the non-local approach. 

The medium modifications of the properties of the nucleon resonance are 
often expressed in terms of changes of its mass and width from their free 
values in a manner dependent on the nuclear density [20, 11]. Our non- 
local approach can accommodate these density dependent parameters for 
the resonance. After implementing the density dependent mass and width 



for the resonance in the s-channel part of the above amplitude, we calculate 
the resulting integral following an approach similar to the one used in ref. [21]. 
We rewrite the integral of the s-channel part in eq(l) as 

(2tt) 4 J d% x t /( i + 1 ^fj i )e ikNy <f>;(y)W(y), (2) 

where we use a plane wave function (see ref. [19]) for the final state of 
the struck nucleon; Xsf denotes its spin state. The function W(y) has the 
following form 

-ip(y-x) 



(2nrW(y) = d*xd*p lM,4^ ^b{x). (3) 

Acting with the operator^ — Mr (p) +i-ip- (p refers to the density of the nu- 
cleus) from the left on (3) and then carrying the integration over momentum, 
the following Dirac-type linear differential equation is obtained 

(j-M R (p)+i^)W(y) = V(v), (4) 



where the source term is 



V(y) = ^Je-^My)- (5) 
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Equation (4) leads to the following second order differential equation for the 
the upper component of W(r). 

[p 2 - a(r)(3(r)]W up (r) = a ■ pV d (r) - (3(r)V up (r), (6) 

where the indices up and d indicate the upper and lower components of 
the functions W(r) and V(r). Note these functions continue to be spin 
dependent. The lower component Wd(r) can be obtained from the upper 
component as, 

(3{r) 
The functions a(r) and /3(r) are 
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(3{r) = E b + w 1 + M R (p)-i 1 ^-, (8) 



where Eb is the energy of the bound nucleon. W up (r) is obtained by solving 
eq.(6) following the same approach as in ref.[15]. 

Substituting the solution of W up (r) into eq.(l) we get, 
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The transition function for the u-channel Z" Mb is calculated in the local 
approximation and is the same as eq.(3) in [19]. The corresponding function 
for the s-channel is calculated in non-local approach described above and is 
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where v\ n (r) is the radial wave function of the 77 meson, kj and k T; are unit 
vectors along the directions of the final nucleon and rj meson, respectively. 
The amplitude (9) is used to calculate the cross section for the exclusive 
reaction (see eq(5) of ref. [19]) . The inclusive cross section is obtained by 
numerical integration of the exclusive cross section over the phase space of 
the final nucleon and summing over all energy levels of the target nucleus. 

In the following section we shall explore whether the available data can 
be used to clarify the issue of possible modifications of the mass and width 
of the Sn resonance. For this purpose we will adopt phenomenological forms 
of these parameters suggested by earlier works [20, 11]. For mass we use a 
density dependent form suggested by the quark-meson coupling model for 
hadrons [20]. 

M(p) = Mfree (l - bjp) (11) 

For width we use a density dependent form given in [11] which includes the 
broadening of the resonance width due to its collisions with hadrons in the 
medium, 

r(p) = r free + c -^, (12) 



3 Results and Discussion 

Throughout the calculations presented here we use scalar and vector poten- 
tials of Woods Saxon form to generate the bound state Dirac wave functions 
of the nucleons. The parameters of these potentials for different targets are 
listed in table 1. The final state interactions of the eta meson are taken into 
account through the use of the optical potential DW1 from ref.[22]. For the 
coupling to the proton we use the parameters of set 3 given in table 1 of 
ref.[23] and in calculating the neutron contributions we use the factor | for 
the ratio of neutron to proton cross sections. 

It is well known that the photoproduction process on the nucleon, in the 
second resonance region, is dominated by the Sn(1535) diagram. We begin 
our discussion by showing that this situation carries over to the inclusive 
reactions on nuclei. The collective contributions from other diagrams are 
small by comparison. This is shown in Fig.l. Figure 1(a) shows the angular 
distribution of the inclusive cross section of the reaction on 12 C at a photon 
incident energy E 7 = 820 MeV. The solid curve gives the Sn contribution 
while the dashed curve includes the additional contributions due to the D13, 
the nucleon pole terms and the meson exchange diagrams. 

The total inclusive cross section on 12 C for photon energies in the range 
0.65 - 1.1 GeV is depicted in Fig.l(b). The dominance of the Sn contributions 
is clearly evident from these comparisons. It is worth mentioning here that 
this rather fortunate circumstance arises not only from the relative smallness 
of the contributions of the other diagrams, but also from strong cancellations 
among these contributions. 

Fig. 2 shows the inclusive angular distributions for C, Al and Cu. For each 
target nucleus, we show the distributions at two photon energies from among 
those measured: 740 MeV and 980 MeV. The data are those of KEK [9], LNS 
[10] and MAINZ (only for 12 C at 740 MeV) [8]. The solid curves represent 
non-local calculations that assume free propagation for the Sn resonance. 
We take the free mass and width of the resonance to be 1535 MeV and 160 
MeV, respectively. In these and subsequent figures these curves are labeled 
'Free'. General trends in the data are reproduced and the agreement is almost 
quantitative in some cases. The main exception is the case of Al at 980 MeV. 
There are some disagreements among the measurements from KEK and LNS 
in the case of Cu and our theoretical calculations show better agreement with 
the latter measurements. 
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The dashed curves in Fig. 2 represent calculations in which, following 
suggestions of ref. [11], the width of the Sn resonance is assumed to undergo 
density dependent broadening in the medium of the form given in eq(12) 
with c = 50 MeV. We note that the broadening leads to reduction in the 
cross sections for all three nuclei. It may be argued that the effect leads to a 
worsening in the accord with the data for C at the lower energy (740 MeV) 
and by contrast to a slight improvement at 980 MeV. For Al the quality 
of agreement with the data remains about the same, whereas for Cu, the 
disagreement between measurements makes it difficult to judge the outcome. 

We now consider the effect of a reduction of mass on the angular distri- 
butions. In the discussion of ref. [20], the authors suggested that the nucleon 
mass in the medium is reduced according to eq.(ll) with a value of the reduc- 
tion parameter b = 0.14. Fig. 3 shows the results of the angular distribution 
for the same targets and photon energies as in Fig.2. The dash-dotted curves 
are obtained when the above value is used for the reduction of the mass of 
the Sn. Compared to calculations with no medium effects (solid curves), we 
note a large decrease in the cross section. This also brings the calculated 
cross section much below the data. This clearly indicates that if the mass of 
the Sn were reduced in the medium, the scale of reduction would be much 
smaller than is indicated by the value of b given above. If the reduction in 
mass is made weaker the effect becomes smaller. For example the dashed 
curves show the calculations for a much smaller value of b (b = 0.014). Here 
the effect on the cross section is relatively small. We note that at the lower 
energy the effect leads to a slight increase in the magnitude of the cross 
section while the effect at the higher energy is a slight decrease. 

The data for the total inclusive cross section (obtained by integrating the 
angular distributions discussed above) are shown in Figs. 4-6. These include 
the data for the three nuclei C, Al and Cu from KEK, Mainz and LNS. Fig. 4 
displays the inclusive data for C. In Fig. 4(a) we show the effects of changes in 
the mass of the Sn on the behavior of the total inclusive cross section in the 
energy range from near threshold to 1.1 GeV. The solid curve represents the 
case of no medium modification of the mass. We note that these calculations 
provide reasonable agreement with the data particularly over and beyond 
the peak, but are slightly lower at the lower energies. The dashed and dash- 
dotted curves represent the two cases where the mass modification parameter 
b of eq.(ll) has the values 0.014 and 0.035, respectively. The latter value, 
though it improves the accord with the data on the lower energy side of the 



peak somewhat, still gives a strong suppression on the higher energy side. 
Calculations with the smaller value of b lie blew the data at energies above 
900 MeV. 

In Fig. 4(b) we present similar comparisons involving changes to the width 
of the resonance. We show two cases for density dependent changes in the 
width of the form given in eq.(12). The dash-dotted curve is for a value of c 
= 25 MeV and the dashed curve is for c = 50 MeV. Both these calculations 
lead to a lowering of the cross section away from the data. 

Figs. 5 shows similar comparisons for the KEK inclusive cross section on 
Al and Fig. 6 compares the calculations with the KEK and LNS data for Cu. 
The qualitative features of the results are similar. Again the data rule out 
any strong reduction (associated with b = 0.14) in the mass of the resonance. 
A reduction in the mass of Sn corresponding to a value of b = 0.014 brings 
the calculations somewhat closer to the data for Al, but leads to a reduced 
cross Section at higher energies for Cu. In the case of the modification of 
width we note that the density dependent width of eq.(12), with c = 50 MeV, 
leads to a suppression of the cross section away from the data. 

We have argued earlier that the dominance of the Sn contributions is 
used to limit the forgoing comparisons and discussions to this contribution. 
In fig 7 we revisit this question. In part (a) we show the results for the 
case of production on C for the Sn. We compare this with the results in 
the lower panel (b) based on calculations in which the small contributions 
from the D13, the nucleon pole and the t-channel vector meson diagrams 
are included. These contributions, being small, are calculated using the 
local approximation (the non-local corrections in this case are expected to 
be negligible). Each of the curves in panel (a) has a counterpart in panel (b): 
The solid curves represent the free calculations, the dashed curves represent 
calculations in which the mass reduction is affected by a value b = 0.014. 
The dash-dotted curves represent calculations in which we make changes in 
both mass and width of the resonance. We assumed the width is reduced to 
150 MeV (from 160 MeV) and in the meantime the mass is reduced as above. 
It is seen from the figures that the inclusion of the additional contributions, 
though it modifies the cross section slightly, does not alter the character of the 
comparison with the data, from those involving only the Sn contributions. 

The preceding comparisons were based on calculations using the optical 
potential DW1 to describe the final state interaction of the 77 meson. The 
calculations of the cross sections are sensitive to the optical potential used. 
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For example we show in Fig. 8 a comparison of the cross sections calculated 
using one of the optical potentials derived by Oset et al. [24]. These calcula- 
tions of the eta optical potential depend on the magnitude of the real part 
of Sn self energy in the medium, and we have used the case for which the 
real part of the Sn self-energy is -50 MeV. The cross sections are seen to lie 
below those calculated using the DW1 potential. It is clear that this type of 
uncertainty should be kept in mind before a final verdict is gleaned form the 
data as to the extent of any medium modifications of the properties of the 
resonance. 



4 conclusion 

A relativistic non-local model for the inclusive photoproduction of 77 mesons 
from complex nuclei is introduced. This model is based on the dominance of 
the contributions due to the formation of the Sn resonance. The inclusion 
of non-local effects is the main improvement over our earlier work which was 
applied to the limited data available at the time. Moreover, the non-local 
approach makes it possible to include mass and width parameters of the 
resonance that may depend on the nuclear density. 

Calculations were first carried out using the free mass and width for 
the resonance (i.e. no medium effect). Comparison of the results of these 
calculations with the data available for C, Al and Cu, reproduce the over 
all shape of the angular distributions as well as the total cross section data 
up to 1.0 GcV. It is worth noting that our results are in somewhat better 
agreement with the recent data of LNS than with the older data. 

The influence of possible medium modifications of the properties of the 
resonance was also studied. Simple density dependent forms for the mass and 
width of the resonance were used to test the effect these have on the calculated 
cross sections and their accord with the data. Comparison with the data 
showed that reductions in the mass of Sn on a scale similar to that suggested 
for the nucleon leads to cross sections that disagree strongly with the data. 
The data can accommodate much weaker changes in the mass and some 
broadening of the width. Although one can claim some scattered instances of 
improvement with these changes, it is not possible to claim that the data show 
clear evidence of medium modifications of resonance parameters. The lack of 
this evidence is due in part to the occasional disagreement between different 



10 



measurements. There is also the inherent weakness in our present knowledge 
of the final state interaction of the eta meson with nuclei (represented here 
by the optical potential), which adds to the uncertainty. In conclusion we 
would estimate that if any medium modifications of the properties of the Sn 
resonance existed, they would likely be rather small. This is in agreement 
with the general conclusions reached by the authors of refs.[7, 10, 11, 12]. 
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Table Caption 

TABLE 1. Strengths, reduced radii and diffuseness for the relativistic scalar 
and vector mean-field potentials, respectively. 



Figure Captions 

FIG. 1. Differential (a) and total (b) inclusive cross section for 12 C( , y,r]) 
reaction. Solid curves: Sn resonance contribution, dashed curves: Sum of 
contributions from Sn resonance, D 13 resonance, nucleon pole and vector 
meson. 

FIG. 2. Differential inclusive cross section for (7,77) reaction on C, Al and 
Cu nuclei. The graphs in the left and right panels represent the results for 
incident photon energies of E 7 = 740 MeV and E 7 = 980 MeV, respectively. 
Solid curves: calculations using the values 1535 MeV and 160 MeV for the 
mass and width of Sn resonance, respectively (labeled Free in this and sub- 
sequent figures), dashed curves: density dependent width of eq(12) with c = 
50 MeV (labeled T(p)). Optical potential DW1 of ref. [22] is used. Data are 
(KEK) from [9], (LNS) from [10] and (MAINZ) from [8]. 

FIG. 3. Same as Fig. 2, but for the effect of mass modifications. Solid curves: 
free calculations, dashed curves: density dependent mass of eq(ll) with b = 
0.014, dash-dotted curves: density dependent mass of eq(ll) with b = 0.14. 

FIG. 4. Total inclusive cross section for the 12 C(7,r]) reaction. Data and 
potentials are from the same references as in Fig.2. Fig.4(a); solid curve: free 
calculations, dashed curve: density dependent mass of eq(ll) with b=0.014, 
dash-dotted curve: density dependent mass of eq(ll) with b=0.035. Fig.4(b); 
solid curve: free calculation, dashed curve: density dependent width of eq(12) 
with c = 50 MeV, dash-dotted curve: density dependent width of eq(12) with 
c = 25 MeV. 

FIG. 5. Same as Fig. 4 but for Al. In Fig. 5(a): solid curve: free calculations, 
dashed curve: density dependent mass of eq(ll) with b=0.014, dash-dotted 
curve: density dependent mass of eq(ll) with b=0.14. In Fig.5(b): solid 
curve: free calculations, dashed curve: density dependent width of eq(12) 
with c = 50 MeV, data are those of ref. [9] . 
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FIG. 6. Same as Fig.5 but for Cu. Curves are labeled same as Fig. 5. Data 
are: (KEK) and (LNS) from refs. [9] and [10], respectively. 

FIG. 7. Same as Fig. 4. Curves in panel (a) are contributions from Sn 
resonance, solid curve: free calculations, dashed curve: free width of 160 MeV 
and density dependent mass of eq(ll) with b=0.014, dash-dotted curve: free 
width of 150 MeV and density dependent mass of eq(ll) with b=0.014. The 
contributions of nucleon poles, D 13 resonance and vector meson diagrams are 
added to the ones for Sn in the calculations presented in panel(b). Curves 
are labeled as in panel(a). 

FIG. 8. Same as the free calculations in Fig.4 (Sn resonance contributions 
only). Solid curve: using optical potential DW1 of ref.[22], dashed curve: 
using optical potential of ref.[24] with the real part of the Sn self-energy set 
to -50 MeV. 
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Nucleus 


MeV 


(fm) 


a,, 
(fm) 


v s 

(MeV) 


(fm) 


(fm) 


12 C 


385.7 


1.056 


0.427 


-470.4 


1.056 


0.447 


^ 7 A1 


354.1 


1.09 


0.450 


-444.5 


1.09 


0.500 


63 Cu 


348.1 


1.149 


0.476 


-424.5 


1.149 


0.506 
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